A TOPOLOGICAL DEGREE COUNTING FOR SOME 
LIOUVILLE SYSTEMS OF MEAN FIELD EQUATIONS 
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Abstract. Let A = (a;j) nxn be an invertible matrix and A" 1 = 
(a ,;, )„ x „ be the inverse of A. In this paper, we consider the general- 
ized Liouville system: 

(0.1) A gUi + a*iPi (jj^— ~ = in M - 



where < hj € G (M) and pj £ R + , and prove that, under the assump- 
tions of (Hi) and (H2) (see Introduction), the Leray-Schauder degree of 
(|0.1[) is equal to 

(- X (M) + 1) ■ ■ ■ (-x(M) + N) 
N\ 

Hp— (pi, ■ ■ ■ , p„) satisfies 

n n 

STYN^Pi < a HPiPj < 8lV ( N + !) Pi- 

i — l l<i,j<n i—l 

Equation (|0.1|) is a natural generalization of the classic Liouville equation 
and is the Euler-Lagrangian equation of Nonlinear function <1> P : 

0/ " 'V.,.// • V„.<; ; • V / /vui-VVlog / foe Ui . 

2 Ja/ l<i,i<n i=l Jju i=l Jm 



The Liouville system (|0.1[) has arisen in many different research areas 
in mathematics and physics. Our counting formulas are the first result 
in degree theory for Liouville systems. 



1. Introduction 

Let (M,g) be a compact Riemann surface with volume 1, hi,...,h n be 
positive C 1 functions on M, pi,..,p n be nonnegative constants. In this 
article we consider the following Liouville system defined on (M,g): 

(1.1) X,,,, ■ }2f>,a,,i k h Je J jdv - 1 ) = 0> i€l:={l,..,n} 
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where dV g is the volume form, A g is the Laplace-Beltrami operator, in local 
coordinates it is of the form: 

When n = 1, equation (jl.ip is the mean field equation of the Liouville type: 

when an = 1- Therefore, the Liouville system (jl.ip is a natural extension 
of the classical Liouville equation, which has profound connection with ge- 
ometry and physics, and has been extensively studied for the past three 
decades. 

If u is a solution of (jl.ip . then after adding a constant, u + c is also a 
solution of (II. ip . Hence, we can always assume u G H 1 (M), where 

fl-^M) = ju G L 2 (M) |V ff it| G L 2 (M) and J udV g = 

For any p = (pi, • • • , p n ), pi > 0, let <£ p be a nonlinear functional defined in 
H 1 ' 71 = H 1 (M) x ■■■ xH\M) by 

2 ijei Jm jei Jm 

where ( 

fl^)nxn is the inverse of A — (^y )nxrc 

It is easy to see that equation 

(jl.ip is the Euler-Lagrangian equation of <P p . 

For a bounded smooth domain $7 in M. 2 , we are also interested in the 
following system of equations: 

(L3) | A "- • >;; : .''^jjt^- ° 

1 Ui\an = 0, i G / 

where /ij are positive C 1 function on £1. 

The Liouville equation (jl.2p or systems (jl.ip and (jl.3p have appeared 
in many different disciplines in mathematics. In conformal geometry, when 
p = 8-7T and M is the sphere § 2 , equation (jl.2p is equivalent to the famous 
Nirenberg problem. For a bounded domain in M 2 and n = 1, (jl.3p can 
be derived from the mean field limit of Euler flows or spherical Onsager 
vortex theory, as studied by Caglioti, Lions, Marchioro and Pulvirenti [3 [8], 
and Kiessling [26\ , Chanillo and Kiessling [9] . In classical gauge field theory, 
equation (jl.ip is closely related to the Chern-Simons-Higgs equation for the 
abelian case, see [BJ [22l [23J, HO]. Various Liouville systems are also used 
to describe models in the theory of Chemotaxis [161 [25] , in the physics of 
charged particle beams [U [191 [27l 128] , in the non-abelian Chern-Simons- 
Higgs theory [20l [24J, HO] and other gauge field models [HI [29] . For recent 
developments of these subjects or related Liouville systems in more general 
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settings, we refer the readers to[Ill2ll3lEllI2llI3llEllI5linilIBll3niEIll32l 
[331 Ell Ell [Ml [Ml El [42] and the references therein. 

For a bounded smooth domain Q in M. 2 , Chipot-Shafrir-Wolansky |17] 
considered equation (|1.3p . where the constant matrix A = (aij) nxn satisfies 
the following condition: 

(HI): A is symmetric, nonnegative, irreducible and invertible. 

Here A is called nonnegative if ajj > for alH, j 6 I = {1, 2, • • • , n}, and is 
called irreducible if there is no subset J of I such that 



In another word, equation (|1.3p can not be written as two de-coupled sub- 
systems. In [17], the authors introduced nonlinear functions Aj(p) of p = 
(pi, ■ ■ ■ ,p n ) defined by 



Among other things, Chipot-Shafrir-Wolansky [T7] proved the following the- 
orem. 

Theorem A. Suppose A satisfies (HI), hi, .., h n are positive C 1 functions 
on Q,, and p = (p%, • • • , p n ) satisfies 



Then equation (|1.3p possesses a solution. 

We note that in [T71 [18] , the authors also proved that the sufficient con- 
dition (jl.4p in Theorem [A] is also a necessary condition for the existence of 
equation (jl.3p when VL is a ball. When n = 1, equation (jl.3p with the param- 
eter p € T is equivalent to the Liouville equation with p < 8n. From various 
known results of Liouville equations, we expect that solutions of equation 
(jl.3p with p r, should have Morse index bigger than 1. Therefore, the 
classical Leray-Schauder degree theory is a suitable tool to be applied for 
studying equation (jl.ip or (|1.3|) when p T. 

To apply the degree theory, we should first prove the a priori bound for 
non-critical parameter p, or equivalently, study the asymptotic behavior of 
bubbling solutions. In |33| . we have proved that near each blow up point, the 
behavior of these bubbling solutions can be controlled well by the standard 
bubble, under the assumption that all the components Uj of solutions have 
to blow up simultaneously i.e., a suitable scaling of Uj should converge to 
some entire solutions of Liouville system: 



djj = for alH G J and j £ I \ J. 




for any non-empty subset J of I = {1, 2, • • • , n}. Let 

r = {{pi, ■■■ ,p n )\pi>0, Aj(p) > for all C J c /}. 



(1.4) 



per. 



(1.5) 



Ar, • V ; , ■ =0 in R 2 

f R2 e idx < +00, i £ I. 
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This blow-up is called type 1. However, there might be situations that at 
some blow-up points, only part of the components of u% (i £ I) blows up, 
but the remaining part does not blow up. This blow-up is called type 2. 
If both type 1 and type 2 occur simultaneously for a sequence of bubbling 
solution, then the set of critical parameters could be very complicated. For 
example n = 2, the critical parameter might be 



{(pi, P2) I Pi = Pi + 8irm, p 2 = p2 + 8irl and (pi, p 2 ) satisfies 
87rfc(pi +P2) = m.l,k G N}, 



where an = 022 = 1 is assumed, N is the set of all natural numbers. Thus, 
the topological degree would be very difficult to compute. In this paper, we 
will prove that this complexity can be avoided if we assume the coefficient 
matrix A to satisfy (HI) and the following (H2) condition: 

(H2): o" < for i G I, a* > for % / j, i,j G I, 
and > for i G /. 

Throughout the paper, we assume that A satisfies both (HI) and (H2). 
For n = 2, A = ^ a |* 022) sa ^ s ^ es C^^) ano - (H2) if and only if aij > 0, 

(0 a\ a 2 N 

01 03 

a 2 a 3 

Then A satisfies (HI) and (H2) if and only if > and a\ + a,j > for 

To state our results, we begin with equation (ll.ip . 

Theorem 1.1. Let A = (aij) nxn satisfy (HI) and (H2) ; and N be a 

nonnegative integer and 

On = {(pi, ■■■ ,p n )\pi>0,i G I and 

8vriV Pi < a vP*Pi < 87r ( N + !) X] 

Suppose hi G C l {M) is positive and K is a compact subset of On- Then 
there exists a constant C such that for any solution u = (u\,--- ,%) of 
(jl.ip with p G K , we have 

\ u i( x )\ — C f or i 6 I an d x ^ M. 

Note that the set Oat is bounded if an > for all i, and Oat is unbounded 
if ajj = for some %. By Theorem II .1^ the critical parameter set for (II. lh is 
the set Tn, where 

T N = < p I MN Y Pi = Y ^iPiPi 

\ iel i,j£l 
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After Theorem 11.11 for p £ Tn for any positive integer N, we can define 
the nonlinear map T p = (T 1 , • • • ,T n ) from H^ n = ^(M) x • • • x H l {M) 
to H x ' n by 

Obviously, T p is compact from H ' n to itself. Then we can define the Leray- 
Schauder degree of equation (jl.ip by 

d p = deg(I -T p ; B R ,0), 

where R is sufficiently large and Br = {u j u € H l,N and \\u\\ < R}. By the 
homotopic invariance and Theorem 1 1.1| d p satisfies the following properties: 

(i) d p is a constant for p £ CV, 

(ii) d p is independent of h = (hi, h>2, • ■ ■ ,h n ). 

The following result is the formula for computing d p . 

Theorem 1.2. d p be the Leray-Schauder degree for U.l\) for p G CV_i ; 
N £ N. Then 

(1.6) d p ={ ! 



where xm is the Euler characteristic of M . 

Since xm = 2 — 2g e where g e is the genus of M, the following existence 
theorem is implied by Theorem 11.21 

Theorem 1.3. (Main Theorem) Let M be a compact Riemann surface with 
genus greater than and hi,..,h n be positive C 1 functions on M. Then 
U.l\) always has a solution for p £ Tjv for any N € N. 

Similarly, for equation (|1.3|) . we have the following result: 

Theorem 1.4. Let (h\,..,h n ) be positive C 1 functions on Cl. Then the 



Leray-Schauder degree d p for (jl.3p is 

dp= [w. ((-X + + N) \ , p G On, N € N 

where x = 1 — 9e, 9e is the number of holes inside O. In particular if $7 is 
not simply connected, (|1.3p always has a solution for p. 

The organization of this paper is as follows. In section 2 we mainly 
address entire solutions of (|1.5j) and some important properties implied by 
the assumption (H2). Then in section 3 we give a detailed description on 
the asymptotic behavior of blowup solutions near a blowup point. In [33] 
the authors have proved that if the system all converges to an entire system 
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of n equations around each blow-up, then all the blow-up solutions converge 
to the same system after scaling. In this section we consider the case of the 
type 2 blow-up, and give a sharper estimate for the bubbling solution near 
the blow-up point. In section 4 by using the sharper estimates in section 3, 
we will prove that type 1 and type 2 blow-up can not occur simultaneously 
in any sequence of bubbling solutions, which leads to the proof of Theorem 
11.11 and Theorem 11.21 in section 5. Finally in section 6 we prove Theorem 
PI 



(2.1) 
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2. On entire solutions 
In this section, we discuss the entire solutions of the Liouville system 
A ^ + Ej 6 7 Oije Ui = 0, M 2 , iel 

J R2 e Ui < oo, iel. 
System (|2.ip is closely related to the following system of equations 
( Avi + //,(— '"••* ' =0 R 2 , iel. 

(2-2) 

I Ji^e^' 6 '^ < oo, iel, 

which was studied initially by Chanillo and Kiessling |10j . and by Chipot- 
Shafrir-Wolansky |17} [18] . Obviously, these two systems are equivalent if 
the coefficient matrix A = (a^) is invertible. In this section, the coefficient 
matrix A is assumed to be symmetric and nonnegative, but not necessarily 
invertible. After a permutation of rows, A can be written as 

/ M \ 



(2.3) 



A 



2 



where each A[ := (ajj)jj e / ; (I = 1, • ■ • , k) is irreducible and / = U/=i^- 
For a positive vector a = {a%, • • • , a n ) (which means each cji is positive), we 
define 

Aj(cr) = <Ji — a ij&i&j f° r an Y C J C /. 

Then Chipot-Shafrir-Wolansky prove the following theorem in |17j : 

Theorem B. (Theorem 1.4 of [17] ): Let A be a nonnegative, symmetric 
matrix that satisfies \2. 3\) . Let a = (01, • ■ ■ ,a n ) be a positive vector such 
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that 

(2.4) A Il (a) = 0, Aj{a)>0, |CJCJ,, l = l,---,k. 

Then there exist v = (v\, ■ ■ ■ , v n ) and [i = (hi, • • • , /u n ) (m > 0) such that 
f2.2\) is satisfied and 

(2.5) °i = ^- [ /<,<-■ 1 i€l. 

Conversely, for an entire solution v to A2.2\) , a defined by \2. 5\) satisfies 



We note that if one submatrix A\ is zero, then I[ consists of only one ele- 
ment because A\ is irreducible. In this case, no positive a satisfies A/ ( (<t) = 
4a = 0. Therefore (|2.2|) has no solution for any positive /i = (fj,i, ■ ■ ■ , fi n ). 

Let Oi (i € I) be positive such that a = (a±, • • • , a n ) satisfies (12. 4j) . Then 
by Theorem B, there is a solution v = (v%, ■ ■ ■ ,v n ) of (|2.2p such that (|2.5p 
holds. In [17], Chipot-Shafrir-Wolansky proved the asymptotic behavior of 
v at oo : 

(2.6) Vi(x) = -<Ji log \x\ + 0(1), iel, 
and 

(2.7) a i3 -(7j > 2, 



due to exp [Ylj&i a ij v j) e ^> 1 (K n ). For any solution t> of (12. 2p . set 

(2.8) Ui = a^Uj + log/xj. 

Then Uj is a solution of the system (|2.ip with 

(2.9) Uj = — nii log |ac| + O(l) for large 
and 

(2.10) mi > 2, 
where 

(2.11) uij =) aijaj. 

Conversely, let u = (ui, ■■■ , u n ) be a solution of (12. ip and set 

(2.12) a i = — e u 'dx is finite, i € I. 

2vr ./ra2 



As did in [17], «j satisfies the asymptotic behavior of (12. 9p where m,j is 
defined by (LTTT]) . By the fact e Ul € L^M 2 ) and by using the Brezis-Merle 
type argument (see Lemma 4.1 in [33]), it can be proved that m £ L% C (R 2 ) 
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and then Uj G C°°(]R 2 ) by further applications of standard elliptic regularity 
theorems. Set Vi to be 

(2.13) Vi ( x ) = — [ log— IfJ—e^Wda;. 

27r J R 2 |x - y| 

Then — ^2j^j a ij v j is a harmonic function in M? and is bounded from 
above by clog|x| for large x. Thus, Ui(x) = Yljei ai i v i^ x ) + Ci ^ or some 
constant q. Clearly, (ui,-- - ,v n ) satisfies f|2.2j> : 

-Avi = e u ' = //, ( >: 

for fa = e c \ Therefore, a = (a±, ■ ■ ■ ,a n ) satisfies (|2.4p . 

We note that in [10] and [17], it has been shown that any component 
Ui of a solution u of (|2.ip with e Ui G L 1 (M 2 ) must be radially symmetric 
with respect to some point in M?, in paticular, if A is irreducible, then m, 
i G I, are radially symmetric with respect to one common point. If A is not 
irreducible, then by (|2.3|) . u«, i € Ii, is symmetric with respect to some 
for Z = 1, 2, • • • , k. By replacing m by 1^(2 + p^), i G 1^, we conclude that 
for any cr satisfying (|2.4p . there exists a radial solution iij of (|2.ip such that 
ai = 2?F /r 2 eUi dx. In summary, Theorem B can be written as the following 
result for system (|2.ip . 

Theorem C: Let A be the same as in Theorem B. Then a solution u = 
(ui, ■ ■ ■ ,u n ) of ()2.ip exists with o~i = j- J R2 e Ui dx if and only if a satisfies 
(|2.4p . Furthermore, after a translation, the solution u is radially symmetric 
with respect to the origin. 

Set 

(2.14) £ := {a; > 0, i G I; A/(cr) = 0; Aj(a) > V0 C J C /.}. 

Let u be a radial solution with a £ £, where a is given by (|2.12p . Then it is 
clear to see that 

Ui{x) = Ui(Sx) + 2 log 6, i G /, 

is also a solution of (12. ip with the same cjj. Thus, without loss of generality, 
we may assume u satisfies itj(0) = an, u n (0) = 0, 1 < i < n - 1, and let 
C IR n_1 be the set of initial values of solutions in (12. ip . Assume (HI) 
holds. In [33J, the authors proved 

(2.15) B is open and is homeomorphic to £ (defined in (I2.14p ). 

In particular £ is an open set in {a; A/(cr) = }. If an > for all i G /, 
then it is not difficult to see B = M ra_1 . In general B might not be equal to 
M n_1 , it is even not known whether B 7^ 0. 

In the following Lemma 12.11 and Theorem 12.11 we show that with (H2), 
all submatrices of A are irreducible, B 7^ and the condition Aj(cr) > is 
automatically satisfied if A/(<r) = 0. 

Lemma 2.1. Let A satisfy (HI) and (H2), thena^ > andm&x(au,ajj) < 
dij for all i j^z j in I. In particular, all submatrices of A are irreducible. 
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Proof of Lemma I2.lt 

Step one: If aj j > for some iq, then aj M > aj j for all j ^ zq. 

Suppose an > 0, let u\ = 1, cr 2 = .. = a n = and rrii = YTj=\ a n a r 
Then rrij = a\j for all j G /. Let m = min{m2, m n }. We want to show 
that m > 0. Indeed, if m = 0, we let 

J = {i G J; rrij = 0.}. 

Clearly 1 J, so for any i G J, <7j = 0, which reads 

= o-j = ^ o y 'm 3 - + ^ a^mj = ^ a^rrij, Vi G J. 
jeJ j£J j<£J 

Since a 13 > for all i ^ j (see (H2)), we have a tJ = for all i G J and j J. 
After a permutation of the rows of A" 1 (therefore the same permutation on 
the columns) A' 1 is of the form: 



Bi 
B 2 



which means A is of the form 





B 2 X 

after a permutation of its rows and columns. This is a contradiction to 
the irreducibility of A. Therefore we have proved m > 0. Next we claim 
m > mi. If this is not true, we have m < mi and we let 

J := {i G /; m; = m }. 

Then by our assumption, J / and 1 J. Thus, for any i G J, the fact 
<7j = yields 



(2.16) = <Ti = ^2a Lj m j + ^ 



> m ^ a l ° ' + m ^ a iJ 
= m^a ij >0, ViGJ. 

Note that we have used (H2) in both inequalities. We see that the second 
line of ()2.16p is a strict inequality unless a*- 7 = for all i G J and j J. 
Thus, (|2.16p yields a 4 - 7 = for i G J and j $ J, a contradiction to the 
irreducibility of A. Step one is established. 

Step two: a-ij > for all i ^ j 

We prove step two by contradiction. Suppose there exist i ^ j such that 
aij = 0. By step one, we have an = ajj = 0. Without loss of generality 
we assume i = 1 and j = 2. We can also assume ais > 0, because the 
invertibility of A implies an > for some i > 3. We can apply a permutation 
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on A to move the positive entry to the third row. Thus the matrix A is of 
the following form after a permutation of rows and columns: 



f 





a 13 








«23 


a 13 


Gt23 


a-33 



\ '■■ '■■ '■■ '■■) 

Let a\ = 1, <T2 = 0, cr 3 = 1 and a 4 = .. = a n = and rrtj = Y^=i a ij a j- 
Clearly mi = 013, 7773 > 013. Let m = min{m2, 1714, .., m n }. We first claim 
m > 0. Suppose this is not the case, let J = {i G I; rrij = 0}. Obviously 
1, 3 J, which implies CTj = for all j G J. Using the same argument as in 
step one, we have a 13 = for all % € J and all j J, thus a contradiction to 
the irreducibility of A. 

Next we claim m > mi. Suppose this is not true, then m < mi. Observe 
that m.3 > m i so m < m3. Let 

J = {i G J; mj = m}. 

Then 1, 3 J, so Cj = for all i G J. Then (12361) yields 

= Oi > m V" a ij ' + 2 (mi - m)a 13 > 0, 

which is a contradiction. Therefore m > m\ is proved. In particular, 7772 = 
023 > mi, which gives 023 > (I13. Since we can switch the first two rows of 
A the same argument gives a 13 > 023, consequently 013 = 023. We use 013 
to represent any nonzero entry on the first row, so we have proved that the 
first two rows are identical, a contradiction to the invertibility of A. Lemma 
12. H is established. □ 

Theorem 2.1. Let A satisfy (HI) and (H2). Suppose a = (ai,..,a n ) has 
positive components and A/(<r) = 0. Then Aj(cj) > for all C J C I. 

Proof of Theorem I2.lt 

First we assume 6 7^ (which will be proved in Lemma |2.2[) and prove 

(2.17) If a € d£ , then 3J C I, such that a { > Vi G J; 

<Tj = Vi ^ J; and a satisfies Aj(cf) = 0. 

Proof of (IgHTj) : 

Let cr fc be a sequence of points in £ that tends to a G d£. Let u fc = 
(ui,..,7i^) be global, radial solutions that correspond to a k . Without loss 
of generality we assume iti(0) = maxj g /7i^(0) = 0. Since e u » < 1, by the 
standard elliptic estimates, there exists a subsequence of u k (still denoted 
by u k ) such that parts of u k converge in C^ C (M. 2 ). There are two cases of 
the convergence of u k to be discussed separately. 

Case one: u k converges to a global solution v = (vi, .., v n ) which satisfies 
the system (|2.ip of n equations. 
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We claim this case can not happen. Indeed, let 



of := — / e u i and <Xi := — / e Vi , i e I. 
2tt J-R2 2ir 



Clearly 



(2.18) 5j := lim of > cr*. 

k— >oo 

Note that a = (oi,..,o n ) and a := {p\, .., a n ) both satisfy Aj = and 
we also have Y^jei a ij a j > 2 for i € / because J R2 < oo. Taking the 
difference on the two equations Aj(cr) = and Aj(ct) = 0, we arrive at 

(2.19) ai i di ~ 2 )^' ~ + a v <7< ~ 2 ^ ~ a ^ = °' 

For each i E I, Y^jei a ij^j — Z~2jei a ij a j > 2. Combining this fact with 
(|2.18p and (|2.19p we have cr, = <7j for all i E I. Thus a E S, which is an 
open subset of the hypersurface Aj = (see (J2.15P ). a contradiction to the 
assumption that a E 88. 

Case two: There exists K C I such that u k converges for j € K and 
Uj(r) —> — oo uniformly in any bounded set of [0, oo) for j K. 

Let / = |i^|, clearly uf must converge, so without loss of generality, we 
assume that the first I components of u k converge to v = (v±, ..,vi) which 
satisfies 



l 

(2.20) - Av l = ^2a ij e^, i = l,..,l, in M 2 

3=1 

and it is easy to show 

Oi := lim of > a, L := — [ e v \ i = l,..,l. 

k— >oo Z7T 

Since J R2 e Vi < oo, by (12. 7jl . we have 

i 

(2.21) ,i, l n J > 2. i=l,..,l. 

3=1 

Let a := (ax, .., oi, 0, .., 0) and cf = (ai,.., cr n ) be the limit of a k . We have 
proved that Ui > ai for 1 < i < I. Although (aij)ixi ma Y n °t be invertible, 
we still have the Pohozaev identity: 

i l 

(2.22) ^2 aijOiOj = 4y~Vj. 

i,j=i i=i 

See Theorem C. Now we claim 

(2.23) di = tJi, i = l,..,l; &i = for i>L 
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Clearly <|2T7j> follows from (g55j) and (12321) . The proof of relies 
heavily on (H2). Set 

n 

mj = a ij a ji i €f J. 

3=1 

We want to prove 

(2.24) rrii > 2 for all i G /. 

Obviously this is true for i = 1, .., I, so we need to prove mj > 2 for i > i. 
To see this, we first observe that m 8 > for all i because a^- > 0, i ^ j, by 
Lemma 12.11 

The proof of rrii > 2 can be obtained by the similar argument of Lemma 
ED Let 

m = min{m/ + i,m/ +2 , ..,m, n }. 

Suppose m < 2, then let J = {i E I; rrii = m} and m* = mmi £ jc{rrii}. 
Same as before, 1, G" J, and Oi = for all i G J. So, (12.16j) yields 

= 0* > m aij + ( m * ~ m ) ^ - °> 
j€l j<£J c 

and then a'-' = for i £ J and j G" J c , which yields a contradiction to the 
irreducibility of A. Hence, rrii > 2 is proved. 

Now we finish the proof of (|2.23p . Certainly, (|2.22p can be written as 
A/((t) = (the last n — I components of a are 0). a also satisfies A(<r) = 0. 
Besides we have di > Oi and m^m, > 2 for all i £ /. Using (|2,19p . we 
obtain <jj = Oi for all i. (I2.23D and (|2.17j) are established. 

Finally, for any a on the surface Aj = with all the components positive, 
we claim that a G £ . Suppose a £. Let oe G £, since Aj(cr) = is 
connected, we can find a path r(i)(0 < t < 1) on Aj = that connects <r 
and <te (r(0) = <T£,r(l) = a). Here we require all the components of T(t) 
(0 < t < 1) be positive. Because a £, there exists to £ [0, 1] such that 
r(to) G d£. But it yields a contradiction to (|2.17p because no component of 
r(to) is zero. Theorem 12.11 is established. □ 

In the proof of Theorem 12.11 we assumed £ 7^ 0, which is established in 
the following lemma. 

Lemma 2.2. £ is not empty. 

Proof of Lemma 12. 2t 

Let £i = ]Cj=i a lJ for i £ I. From (H2) we know £j > for all i G /. 
First consider the case that all £j are positive: 

Case one: £j > for all i E I 
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Using the properties of inverse matrices 

n n n n 

(2.25) Yl = E ai i E ajk = E <^ = L 

3=1 j=l k=l k=l 

Let U{ = 4£j for i & I. Then by the assumption in this case Oi > for all 
i 6 I and direct computation shows 

(2.26) Aj(<t) = 4 <7i - ^ Oij-o-iOj- 

= i6(2 & - E = 16 E " E a ^i) = °- 

iS-T ijGi" igl j=l 

For any nonempty J C J, without loss of generality J = 1, for some 
I < n, easy to see 



i=l i,j=l 

I I 

= 16 (E & ~ E "'.^'S/ ) 
i=l i,j'=l 

II In 

= i6j>(i - E a ^) = 16 E( E 

i=i j=i i=i j=i+i 

Clearly Aj > 0. Since by Lemma l2.1| q,j > for 1 < i < I and I + 1 < j < n 
we have Aj > 0. Therefore a € £ and £ 7^ 0. 

Case Two: There exists £j = 

First we observe that it is not possible to have all £j = because otherwise 
adding all the rows of A~ l to the first row would make all the entries of the 
first row 0, a contradiction to the invertibility of A -1 . Without loss of 
generality we assume = ... = £ n = and £1, > 0. Let J = {1,..,/} 
then we claim 

(2.27) Aj(cr) = and A Jx > V0 C J x C J. 

From (|2.26p we see easily Aj = 0. For C Jx C J, without loss of gen- 
erality we assume J\ = with l\ < I. Similar to case one (using 
Ylj=i a ij£,j = 1 an d &i = 4£j for all i = 1, .., Z) 

A ^ = 16 E( E tt ^)&- 

i=l j=h+l 

Thus Ajj > is an immediate consequence of Lemma 12.11 ()2.27p is estab- 
lished. 
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Let A = (ajj)jjgj. Although A may not be invertible, Theorem B can 
still be applied to conclude that there exists a radially symmetric solution 
u(r) = («i(r), -,ui(r)) to 



(2.28) 



Au, t + Yl j= i oye^ = 0, M 2 , i = 1, .., I, 



2^ Ir 2 e " ! = a i = 4£i> * = 1) -J ^ 
Since e ni G L 1 (M 2 ), we have 

/ 

(2.29) X) a y' cr J > 2 ' i = 1 > 

It has been discussed in the proof of Theorem 12.11 that (H2) implies 

I 

(2.30) a n G i > 2 ' Vi e J - 
which leads to 

(2.31) -!- I Vdj/'>2 + ^ iG7 

for some 5 > and R > 1. Now we construct a sequence of functions 
-u e = (tif , as follows: 

f Kr^ + ^O'^ + E^i^^^^O, 0<r<oo, iel, 
uj(0)= Ui (0), i = l,...,l, 

k u|(0) = loge, i = / + 1, ..,n. 

By standard ODE existence theory, solution u € is well defined for all r > 0. 
Since ii|(r) are decreasing functions, it is easy to see that as e tends to 0, 
u\ {% > I) tends to — oo over [0, R] for any fixed R > 0. Therefore, the first 
I components of u e converge uniformly to the corresponding components 
of u over any fixed [0, i?]. In particular for the R in (|2.3ip we have, for e 
sufficiently small 

(2.32) -L f V aije"! > [ V aye tt S > 2 + 5/2, Vi € J, 

27F J Br 27r ./Bh ^ 



which implies 



(«f)'(r) = --/ > ^a^e^^sds < r > 12. 



J ° .7=1 



Thus 

*4( r ) < -( 2 + V 2 )!og r + 0(1), for r > 1? 
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where 0(1) is a constant independent of r but may depend on e. Hence for 
e > small, L 2 & u% < oo for all i e I. Thus £ ^ 0. Case two and Lemma 
J\ are established. □ 



Remark 2.1. For the 2x2 case, let A = ( ° ^ j be an non-negative and 

irreducible matrix (which means c > 0). It can be proved that the conclusion 
of Theorem \2.1\ holds for A if and only if max(a, b) < 2c. However, A 
satisfies (HI ) and (H2 ) if and only if max(a, b) < c and c 2 ^ ab. 

3. The case of partial blowup 

In this section we consider the following case: Let u k = {u k , ..,?%} be a 
sequence of solutions of 

n 

(3.1) -Aii* = Y j a ij h){x)e uk i , B u iel 

j'=i 

where h k = (h k , .., h k ) is a sequence of positive functions with uniformly 
bounded C 1 norm: 

(3.2) 1 < h k (x) < C, \Vh k \(x)<C, Vi G Bi, ViG/. 
Suppose that u k blows up only at 0: 

(3.3) M k = max v%(x) -> oo, max?/* < C(A) VJf CC B x \ {()}. 

i£l,x£B\ i£l 



Let 

(3.4) o* 



— / ^e tt i; m^J^ayaf, iel. 



Without loss of generality, we assume u k (0) = M k , and set 

(3.5) v k (y) = u k {5 k y) + 21og,5 fe , where 5 k = e~^\ 

By elliptic estimates, we can show that there is J C I such that {u^}j e j 
converges, in C 2 oc (M. 2 ), to {t>i}igj. We may assume J = {1, 2, • • • , I}. Then 
v = (vi, ■ ■ ■ ,vi) solves the following subsystem: 

I 

(3.6) - Avi = "■P( ! '''' , i = !,■■■ J, in I 2 , 

where Hi = lim^oo h\ (0). In this case, we have v k (x) — > — oo for i J as 
fc — >■ +oo in any compact set of M 2 . 

We also assume u k to have bounded oscillation on dB\ : 

(3.7) \u k (x)-u k (y)\ <C, Vx,yedB u iel 
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and uniformly bounded energy in B\. 

(3.8) f h\e u i <C, i G I. 

Jb 1 

When I = n, it was proved by the authors {33] that there exists an entire 
solution U k = (Uf, ■ ■ ■ , E7*) of 

Au i +Y1 "'f 11 ! 1 ' ' = in r2 

such that 

\u k (x) - Uf (x)\ < C for \x\ < - and i G I. 
An immediate consequence is the following estimate: 

(3.9) u k {x) = -^-^M k + 0(1) for \x\ = L 



In this section, we want to extend the estimate (|3.9p to the case when there 
are only I components of v k , I < n, which converge to an entire solution of 

dm 



Proposition 3.1. Let h k = (h k ,--- ,h k n ) satisfy jlS.ty) . and let u k be so- 
lutions of Q£Tj} such that f&3|) . (31$ . (3?T\j and hold. Let ct; = 
limfe_ >+QO o\ and mi = linifc_>.-|- QO rnf , where a k and m\ are given by (|3.4|) . 
Then 



(2) cri, .., o { > 0, a l+ i = ... = o n = 0. 

(3) mi > 2 Vie I. 

(4) For 1 < i < I 

(3.10) + ^^M k \ < C(e), Vi 6 Si \ B e , 

(5) For / + 1 < i < n 

(3.11) \u k {x) + ^^M fc + (M k - u k (0))\ < C(e), V* G #i \ S e . 

Proof of Proposition I3.lt 

This proof is divided into a few steps. 
Step one: rrij > 2 for all i & L. 
Let 

0-fo = ^ / Hi(0)e^, i = l,2,..,Z. 
Then Theorem C and (|2.7p imply 

Ei,j=l a ij°~ivO~j v = 4 Ei=l ^iv 



QiijOiO'i . 



i,j=l 



(3.12) 



Zjj=l a ij°~jv > 2, 
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Obviously, Oi > 0{ v for i = 1, .., I. Let <Ji v = for i = I + 1, ..,n and 



3=1 



We claim that > 2 for all i = 1, ..,n. For i < I, this is already known 
by (|3.12p . The proof of rrii v > 2 for i > i is the same as (I2.24D in the proof 
of Theorem 12.11 Since o~i > (?i V , we also have mj > 2. 

Step two: n^(x) — > — oo over any compact subset of B\ \ {0} for all 

i = 1,2, .., n. 

We first show that tends to — oo on dB\. Let G(x,y) be the Green's 
function with respect to the Dirichlet condition on dB\, the Green's repre- 
sentation formula gives 

Ui(x)> J G(x,r?)^^a i j/i^(r/)e n J^ ) ^dr? + rmnnf. 

If mingsi u\ > —C for some C > 0, we use 
z 

aij/tj e u J — > 2-7rmj5o in measure 

3=1 

where <5q is the Dirac mass at 0, and 



G(x, rj) > ——— log |x — r)\ — C\ for \x\ < - 
2tt 2 



to obtain 



lim e u i^ > C2\x\ 2 ei for < \x\ < - and for some t\ > 0, 

k— >oo 2 

where C\ and Ci are two positive constants. This is a contradiction to 
f B e u * < C. Therefore we have proved ming^ — > — oo. Since is 
bounded from above over any compact subset of B\ \ {0}, standard elliptic 
estimate implies that — > — oo uniformly in any compact subset of -E>i\{0}. 

Step three: 07+1 = ... = a n = 0. 

Let wf = YljeJ a ^ u j- Then wf satisfies 

We can apply the Pohozaev identity to the above equation, i.e., multiplying 
x-V(Za 

k — > +oo, 



x ■ ^(Y^aijWj) and applying integration by parts, we obtain by passing 



i,j=l i=l 
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On the other hand, a v = {a\ v , ai v , 0, .., 0) also satisfies the equality above. 
Besides, we also have 

Ci > <?iv, i = l,..,l, Oi > 0, i = l + l, .., n. 

By taking the difference of the above equations satisfied by a and a v , we 
arrive at the same as ()2. 19|) : 

jei" \iG-f / j'Gi" iG-f 

Recall mj = ^ OyOj > 2 and = a,ijdj v > 2. Then the above identity 
leads to a = a v , i.e., <7j = c% for i < / and cij = for i > /. 

Step four: u\(x) + 2 log \x\ < C 

By scaling, it is equivalent to proving 

(3.13) v>°(y) + 21og\y\<C, \y\<8-\ 

Since only {vf,..,vf} converges to {vi,..,v{\ in C 2 oc (R 2 ), it implies that 
vf +1 , ...,t>^ tend to — oo in any compact subset of M 2 . Suppose (I3.13P does 
not hold, there exists y k — > oo such that, along a subsequence 

(3.14) vf (y k ) + 21og \y k \ = max (vf(y) + 21og \y\) -> +oo. 

y|<<5fc 1 ,*eJ 

It is easy to see that \y k \8k —> as is bounded above in any compact 
subset of B\ \ {0}. Consider y € B(y k , \y k \/2), from 

vf(y) + 2 log |z/| < ujd/fc) + 2 log i € / 

we have 

(3.15) vf(y) < v k l0 {y k ) + 21og2 for y € B(y fc , |y fc |/2) and i G /. 
On one hand, the fact — > oo implies 

(3.16) lim ±- / h^(S k -)e^ > a iv , i = l,2,..,l. 

fc^oo 27T 7B(0,|y fe |/4) 

On the other hand, we can set 

wHv) = v i(yk + r k y) -Vi Q (y k ), for |y| < i £ I, 

where r fc = e ~2%te"=) anc i R k = r k \y k \/2 — >■ oo as k — > oo by (|3.14[) . By 
(|3.15p . u>f(y) is bounded from above by 2 log 2, and w^ o (0) = 0. From 
the standard estimate of elliptic equations, there is J C / such that wf(y) 
converges in Cf oc (M. 2 ) for i € J and wf(x) —t — oo for i J in any compact 
set of M 2 . Clearly, io G J. In particular, 



e Wi o > 5 > for some 8 > 0. 



|2/|<-Rfc 
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Thus, together with (13.161) we have 

hi e Ui ° > a k)V + 5 > a iQV , 

Bi 

a contradiction to step three. Hence, step four is established. 
Step five: Estimate of the decay of vf. 

First we choose R » 1 such that 

(3.17) -L f Yl aatfihyV^dy > 2 + 5, i G I. 

Z7r J B R/2 j=1 

for some 5 > 0. Note that (|3.17p obviously holds for 1 < i < I because of of 
the convergence from (vf, ..vf) to (vi, ..,vi). Then by step one, it also holds 
for i > I. In this step we study the behavior of vf for \y\ > R (i 6 I). For 
each R < \y\ < 5^/2, let 

vf(z) = vf(\y\z) + 2log\y\, \<\A< 2- 

By step four vf < C\ for some C\ over B 2 \Bi/ 2 . Consider the equation for 
vf in B 2 \B 1/2 : 



-Avf(z) = ]T aijh^5 k \y\z)e^ B 2 \ B 

Let ff satisfy 



1 . 

2 



-Atf(z) = E"=i <Hjht(5 k \y\z)e v s B 2 \ B 1/2 , 

f?(z)=0, on dB 1/2 UOB 2 . 
Clearly ff > in _E?2 \ -B^- As a result of the upper bound of vf, we have 
(3.18) 0<ff(z)<C 2 , zeB 2 \B 1/2 . 

Obviously, vf — C\ — is a non-positive harmonic function. Hence the 
Harnack inequality holds: 

-mm (vf -d-ff)<C (-max (vf -C x -f* 



Equivalently, 



^t. 1 . 4 „ 
maxv, < —rn.va.Vi + C 2 . 



dBj. 1 C dB x 
Going back to vf, we have 

(3.19) max^(y) < - mmvf + (-2 + -) logr + C 2 , 

oB r G oB r (_j 

for R < r < 57 1 and i = 1, .., n. 
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Let vf(r) = fg B v k be the average of v k on dB r . 



13 r 



By (|3.17j) we have 



2 + 5 



(v?Y(r)< '—, r>R 



for some 5 > 0. So for r > R, 

v k {r) < -(2 + 6)logr + tf(R). 
By (|3.19p . we have, for i € / 

(3-20) ^(y) < I^(r) + (-2 + |) log |y|+C 2 , 

< -(2+ A)log|y| + C 2 + ^(i?). 
Note that for i > I, vf(R) -+ -oo. Thus 

(3.21) ^(y)<_ jRfc _( 2 + A)iog|y|, Vt>i 

for some Rk — > oo as k — > oo. 

Step Six: The proof of G3A0\l and (l3nT]> 

Let 2/1,2/2 be two points on the same circle centered at 0: \yi \ = | z/2 1 • Then 
the Green's representation formula gives 

(3.22) v k ( yi )-v k (y 2 ) 

n 

(G( yi ,r,) - G(y 2 ,i 1 )){Y j a lj h k j {5 k i 1 )e v ")dr 1 

+h k ( yi )-h k (y 2 ) 

where h k = (h k , ■■,h k l ) is harmonic defined by h k = v k on dB(0, 5^ r ). Clearly 

\h k (y)-h k (y')\ = 0(l), V^'eBiO,^) 
because v k has bounded oscillation on dB(0, 5^ 1 ). Next we claim 

(3.23) \v k ( yi )-v k (y 2 )\<C, Vr G (O,^ 1 ), V| yi | = \y 2 \ = r, tel. 

We only need to verify (|3.23p for r < 5^ /2, as the case for r > \5^ is 
obvious. To evaluate the expression for v k (y\) — v k (y 2 ), we first have 

GM ~ Gfar,) = - log |^ + ^ log ]r -^ | 

where yi is the conjugate of j/i when it is considered as a complex number. 
To prove (13.23R we decompose ^(O,^ 1 ) into the following four sets: 

Ei := {v, \v\ < r/2}, E 2 := {rj; \q - Vl \ < \ v - y 2 \, T - < \rj\ < 2r}, 
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E 3 :={??; 1*7-3/21 < \v ~ Vi\, \ < M < 2r}, £4 = 5(0, x ) \ (uf =1 #i). 
By (12301) 

e <*(v)<(l + |y|)-3-* jel, \y\<5^ 

where 5' > 0. Applying this decaying estimate, we can prove the integral of 
the right hand side of (|3.22[) over each Ej is bounded. We omit the proof 
because it is a standard computation. Thus (|3.23p is proved. Once (|3.23p is 
established, we have 

(3.24) v k {y)=v k (r) + 0{l), iel, \y\ = r. 

For vf(r) we have 



1 f n 

(3.25) (v?)'(r) = -— Y, a ^{5 ky )e^dy. 

3=1 

For r > R, the decay rates of v\ in (|3.20|) and (|3.2ip imply 

(M6) h L % - h ( C - L^ B , ) 

= Hl + (r-'- s ). 
r 

Using (13361) in (13351) we have 

(3.27) v k {r) = v k {0) - m k log(l + r) + O(l), < r < S' 1 . 

Clearly (13TT0D and (|3TTD follow from (1331]) and (13371) . This completes 
the proof of Proposition 13.11 □ 

Remark 3.1. It is easy to see from (13.6D that the sequence u k (0) — is 
bounded for 1 < i < I and tends to —oo for I + 1 < % < n. 

4. The strong interaction between bubbles 

In this section, we suppose u k has two blow-up points p\ and p2- By 
Proposition 13.11 at each blow-up point pi, u k after scaling will converges to 
an entire solution of a subsystem (|3.6j) . The following question naturally 
arises: 

Are these two entire solutions the same? 
The following theorem will answer this question affirmatively 

Proposition 4.1. Let Qq be an open and bounded set with smooth boundary, 
Pi,P2 £ be two distinct points. Suppose u k = (u k ,..,u k ) satisfies $3. 
( 3. 7| ) and $3.8\) on and and h k satisfies A3. 2t) over Qq a s well. Let pi,P2 
be the only two blow-up points of u k on Q: 

^%tk ~ * Pti H G I) such that u k (xtk)^oo, t = 1,2. 



maxuf < C(K), VK CC fi \ {pi,P2}, 
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Then for 5 < \ \p\ — P2 \ 

lim / h^e^dx = lim / h^e^dx, i £ I. 

k ^°°JB( Pl ,S) k ^°°JB(p2,S) 

Proof of Proposition I4.lt Let 



of = — / hfe"* dx, a* = — / hfe * dx, 

27r JB( Pl ,S) 2tt Jb{ P2 ,S) 

for 5 > small and i 6 J, Also we let 



i=l 3=1 

We use cjj,mj,CTj and rrtj to denote the limit of af,mf,af and mf, respec- 
tively. Let 

Mfc = maxM^(x), 3; £ 
iei 

and Mfc is attained by some component of u k at p\ k which tends to p\. 
M k an d P2k can be defined similarly. By comparing the value of u\ over 
f\) \ (B(pi, 6) U B(j>2, 6)), using Proposition 13. II we have 

(4.1) + (M fe - u?(pi fc )) = + C^k " «i (P2fc)) + O(l). 

Here we remind the reader that, for example around p±, if the first I com- 
ponents of u k converge to a system of I equations after scaling, then M k — 
u i(Pik) are uniformly bounded for 1 < i < I. In this case, M k — u k (p\ k ) can 
be combined with the O(l) term. For i > I, M k — u k (pi k ) tends to +00. The 
right hand side of (I4.1D can also be understood this way. For each i G /, if 

M k -u k (p lk ) >M k -u k (p 2k ) 

we let 

l k = (M k - u k {p lk )) - (M k - u k (P2k)), I? = 0. 
On the other hand if 

M k -u k (p lk ) <M k -u k (p 2k ) 

we let 

l k = 0, l k = (M k - u k (p 2k )) - (M k - u k {p lk )). 

Set 

/i :={»€/; lim A- > 0. }, 7 2 := {* G I; lim A > 0. } 
and ^3 be the compliment of Ji U J2. From this definition we see easily that 

iini 2 = 0. 

We claim that /1 is empty. We prove this by contradiction. Suppose I\ 
is not empty, then we consider the following two cases: I2 is not empty or 
I2 is empty. 

Case one: I2 ^ 
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Let 



M ft - ft 

A = lim Si = lim Si = lim 



k^oo Mfc k— >oo Aifc fc— >oo Mf. 

We claim that all these limits exist. Indeed, using the definitions of l\ and 
if, (|4.ip can be written as 

™?~2 M fc _ m k i -2 If 

2 M k M k 2 M fc 1 



Take i £ /j, the RHS tends to (r/ij — 2)/2, which implies that along a 

2 



subsequence, the two terms on the LHS are m \ 2 A and 6%. On the other 



hand, take j € /2> the LHS tends to — A, then the RHS has to tend to 



m J 2 2 + (5j. Now (|4.ip can be written as 



(4.2) A — ^ h <5i = — — h (5j, tel. 

From the definition of 5i and Si, we observe that for each i € I\, Si = and 
for z € I 2 , <5, = 0. By (|3. 10[) and (|3.1ip of Proposition 13.11 we have 

(4.3) ai = Q,ieIi; ai = 0,iel 2 . 

Since Si = for i tfL I\ and <7j = for i G /i, we have <7j#j = for all i G J. 
Similarly, cfjJj = for all i £ /. 

Without loss of generality we assume 

/l = {1 < i < *o}- 
For each i E I 2 , the fact crj = yields 

(4.4) = Oi = ^2 a lj fhj = ^2 alj rnj + ^ a ij fhj. 

j£i j£h j02 

We observe that the last term is positive, because fhj > 2 and there exists 
a 13 > for some i £ I 2 and j I 2 . Multiplying Si to the last term and 
taking the summation on i for all i in I 2 , we have 

(4.5) j2(J2 aijfh ^ >0 - 

Combining (|4.4p and (|4.5p . we have 

a lj rhi5j < 0. 

Trivially, there exists i £ I 2 such that 

(4.6) a%- < 0. 

ie/ 2 
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Multiplying a 13 to both sides of (14.2j) (with i replaced by j) and taking the 
summation on j, it leads to 

(4.7) y, ^ (^) A +E a ^ = E & +E^- 

Using the definition of as well as 5% = for i /j , we can write the LHS 
of flUD as 

Since the first term and the third term are both nonnegative, the LHS is no 
less than the second term. Similarly the RHS can be written as 

Note that we have used 5i = for i li- The first term of the above is 
(because i € I2) and the last term is negative (because of (|4.6p ). Therefore 
the RHS is strictly less than the second term. Putting the estimates on both 
sides together we have 

Since ^2j £ ja' 13 > ((H2)) we conclude A > 1. On the other hand, by 
exchanging I\ and I2 in the above argument, we obtain A < 1. Thus we 
have ruled out the first case. 

Case two: I2 = 

One immediately has 5i = for all i 6 I. Hence, the limits A = 
lirrifc-j.oo Mk/Mk and 5k = lim^oo l^/Mk both exist and (|4.2I) holds with 
5^ = 0. Here we recall that both a = (crx, ••,cr n ) and a = (a\, ..,cT n ) satisfy 

which can be written as 

(4-8) £ a-(^)(^) = £ 

We further remark that 

(4.9) a ' lj > 

because for each i, J2jei a ^ — and A~ x is non-singular. To prove our 
result, we need another fact: 

(4.10) All the eigenvalues of F are nonpositive, 
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where 

F = (« y )ioxi hJ e h = {l,.-,*o}- 
Indeed, let fi be the largest eigenvalue of F and rj = (rji, .., r)i ) be an eigen- 
vector corresponding to \x. Here r/ is the vector that attains 

maxv T Fv, v T v = 1. 



veit 



Since a %3 > for all i ^ j, we can choose rji > for all i £ I±. For each 

» e Ji, 

= o-j = V] a^mj + V] a y 'mj. 
jeh 30i 

Plainly by (H2) 

^ a^rrij < 0, ie/i. 

Multiplying 77, on both sides and taking the summation on i, then we have 

> CL^Vimj = Yl t JL Vj rn j- 

i,jeh jeh 

Since each rji > (one of them is strictly positive) and > for i £ I±, we 
have ^ < 0, and it proves (|4.10p . Now we go back to our proof to rule out 
case two. 

Since Si = in (|4.2p . the Pohozaev identity (|4.8p for a can be written as 
£ a«(^A + *)(^A + «,) = £ a* 



Expanding the LHS of the above and using (|4.8p again for Cj, we obtain 

(4.11) A 2 ^ « + 2A aij C^Y^)5j + E aij5i5 i = E QlJ - 

i,jei i,jel ~ ije/i i,j 

The third term of LHS is nonpositive by (|4.10p . The second term of the 
LHS can be written as 

jei iei jei iei 

because crjSj = for all j € I. Thus we conclude from (|4.1ip that A > 1. 
On the other hand from = (i 6 Ii), argued as (|4.6() we obtain an index 
i £ Ii such that 

(4.12) ^ a% < 0. 

3 eh 

Then as we did for (|4.7p . we obtain 

je/ jei jei 
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Following the same calculation as before we obtain 

0>^>% = f + (A-1)J>& 

which forces A to be less than 1. Thus, we have obtained a contradiction to 
A > 1. Case two is also ruled out. Thus we have proved that I\ has to be 
empty. Using exactly the same argument we also have Ii = 0. 

Since l\ = 1% = 0, fj4.2j) becomes 

uij - 2 _ fhj-2 . 
2 ~ 2 ' * 

Using (|4.8p for both (mi,..,m„) and (mi, .., fh n ) we have A = 1. Conse- 
quently <7j = <jj for all i E I. Proposition 14. 1 1 is established. □ 

5. Proof of Theorem 11.11 and Theorem 11.21 
Let u = (ui, ..,u n ) be a solution of (jl.ip . Define 

(5.1) ^ = m - log / hie^dVg. 

Jm 

Clearly v = (yi,..,v n ) satisfies 

(5.2) f hie Vi dV g = 1 

Jm 

and 

n 

(5.3) A gVl + f>.j"<j {ll j<" - 1) = 0, i € J. 

i=i 

To prove the a priori bound for u, we only need to establish 

(5.4) \vi(x)\<C, VxGM, iel. 
Indeed, once we have ()5.4|) . for u we have 

(5.5) log f hie Ui -C < Ui (x)< log f h ie Ui +C, Vx G M. 

Jm Jm 

Since u €i? 1 (M), there exists xo such that u(xq) = 0. Using this in (|5.5p 
we have 

(5.6) - C < log f Ke u ^ < C. 

Jm 



With (|5.ip and (|5.6p we see that a bound for u can be obtained from the 
bound for v. To prove (|5.4p we only need to prove 

(5.7) Vi<C, iel. 

because a lower bound for Vi can be obtained easily from the upper bound 
in (|5.7p by standard elliptic estimate. So we only need to establish (|5.7jl . 

We prove (|5.7[) by contradiction. Suppose there are solutions v k to (|5,3D 
such that maxM,ie/ vf(x) — >■ +oo. We consider the following two cases. 
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Case one: p\ — > p{ > as k — >■ +00 for all i £ I. 
The equation for i> fc is 

n 

(5.8) A s vf + Pjaij(hje v J - 1) = 0, i € J. 

i=i 

In [33] the authors prove a Brezis- Merle type lemma (Lemma 4.1) which 
guarantees that there exists a positive constant eo > such that if 



k 

e" i dx < €q for all i G I, 



B(p,r ) 



then 



(5.9) vf(x)<C, x€ B(p,-£). 

Thus, v k blows up only at a finite set {p\, ■ ■ ■ ,pn}- Since vf(x) is uni- 
formly bounded from above in any compact set of M \ {pi, • • ■ ,Pn}, by 
([5T5]) . v\ converges to Ya=\ m,i(pi)G(x,pi) in C™ C {M \ {pi, ■ ■ ■ ,pn}), where 



(5.10) 



k 

°M) = lim fc^oo h Ib( Pi ,5o) P k j h i eVj dV 9 



for some 5 > such that B(pi,25 ) n B{p Vl 25o) = 0, I ^ /'. Here, G(x,p) 
is the Green function: 

J -A g ^G(x,p) = S p - 1, 
\ / M G(x,p)^(x) = 0. 

To apply Proposition ^. H we rewrite ()5.8p in local coordinates. For p G M , 
let y = (y 1 , y 2 ) be the isothermal coordinates near p such that y p (p) = (0, 0) 
and y p depends smoothly on p. In this coordinate, ds 2 has the form 

e ^ P ) [ (dy l)2 + (dy 2 )2] f 

where V0(O) = 0, </>(0) = 0. Also near p we have 

Ay p (j) = —2Ke^, where K is the Gauss curvature. 

When there is no ambiguity, we write y = y p for simplicity. In this local 
coordinate, (|5.8p is of the form: 



(5.11) - Avf = e*Y^ a ijPj {hj^i - l) in B(0, <5 ), i G /. 

i=i 

Let ff be defined as 

n 

-Aft = -e+ 2 Pj<kj in B (°> 6 o), i € I, 
and /f (0) = |V/f (0)| = 0. Let 5 t * = - and 
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Then (15. lit) becomes 

n 

(5.12) - Avf = ",,//}<'• in B(0, 5 ). 

i=i 

Here, we observe that f B{0jSo) Hfe^dx = f B{0)5o) P^h ie v "dV g . 
Since v\ converges in M \ U*=i B(Pji 2^o)j we have 

AT 

(5.13) |5ftz)-wfty)| <C, Vx,yeM\\jB( Pj ,26o), i € I. 

i=i 

By (|3.1U|) of Proposition we also have 

(5.14) / hid 3 * dV g Q, i€l, 
and by Proposition 14.11 

(5.15) lim / fthie^dVg = lim / p k h ie v Uv g 

k~>o° Jb(p,,6o) k ^°° J B( Pm ,S ) 

for i E I and for any pair of integers I, m between 1 and N. (|5.14p combined 
with ([57T5]) yields for i e / and j € {1, 2, • • • , iV}, 

a, = lim ^ / H^dx = lim J- / p\h^ dV g = 

k^oo Z7T Jb(0,5 ) k ^°° l7T Jb( Pj ,8 ) 27riV 

On the other hand, (a±, ■ ■ ■ ,a n ) satisfies the Pohozaev identity: 

(5.16) 4^<Tj = dijOiOy 

iei i,jei 

Consequently, 

n n 

sttn p,j = aijPiPj. 

Thus, a contradiction to the assumption of the theorem. 
Case two: lim^oo p\ = pi > 0, % = 1, .., Z, lim^oo = for % > I. 
Let M k = maxju^, and M k = {ff +1 , ■•, ^nl- We first show 

(5.17) M k -M k <C 
by contradiction. Suppose M k — M k — > oo, let 

Vf(y)=v?{e~y+ Pk )-M k 

where p k is where M k is attained: vf Q (p k ) = M k . Clearly i$ > I. Thanks 
to the fact that — > — oo for i < I and p\ — > for i > I, converges in 

cL(K 2 )to 

f -AV i0 = 0, M 2 , 
\ 7< (0) = 0, V i0 <0. 
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Clearly Vi = 0, f B e 'o can be arbitrarily large if R is large, this is a 
contradiction to (|5.2p . (I5.17P is proved. 

We use the same notations as in case one. Let pi, ..,pn be blowup points 
for vf. Then around each blowup point, say p±, the equation for v k can 
be written in local coordinates as (j5.12[) with v\ and Hf defined the same 
as in case one. Without loss of generality we assume p\ > for all k and 
I + 1 < i < L and p\ = for all k and all L + 1 < i < n. Then we observe 
from the definition of that — s>0for/ + l<i<L and = for 
i > L. To reduce case two to case one, we need to adjust the terms involved 
with these vanishing H^s. To do this we set /• as 



-A// 



Z^j=L+l a ij e 3 



-M, 



fk = on dB(0,5). 



Since max^ — M/% is uniformly bounded for all i, we have 

\\f?\\d<C 
for some C independent of k. Now we define 



M k +t 



vf + log p\ + 

yk n/T. i ik 



I + 1 < i < L, 
L + 1 < i < n. 



and 



H 



Hh 



J i 



l<i<L 



Hfre-tf = e <P+f?-f? hi, l + l<i<L, 

Pi 



e A , L + 1 < i < n. 
Easy to see there exists c > independent of k such that 
1 



< A? < c, 



|VHf|<c, B(0,S). 



On the other hand satisfies 



% 

-At>* = , B(0, 5), iel. 

Easy to observe that maxu^ — — > — oo for i = Z + 1, n. Therefore case 
two is reduced to case one, which gives 

ViiPt) = °i{p m ) Vi, m € {1, .., N}, VI < i < I. 

Note that &i(pt) = for all z > Z and all t because p\ — > for i > I. Then 
as in case one we obtain a contradiction. Theorem 11.11 is established. □ 
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Theorem 11.21 will be discussed in two cases. 
Case 1. One of an is positive. 

We may suppose an > 0. Thanks to Theorem 11.11 the Leray-Schauder 
degree of (jl.ip for p S On-i is equal to the degree for the following specific 
system corresponding to (pi, 0, .., 0): 



(5.18) 



f M he^dV g 



0, 



AgUj 



J M h ie ^dV g 



1=0 for j > 2. 



where p\ satisfies 
(5.19) 



8tt(7V - 1) < anpx < 8ttN. 



Easy to see (pi,0, ..,0) € On-i- Using Theorem 1.2 of [13], we obtain the 
degree counting formulas (|1.6p in this case. 



Case 2. an = for all i € /. 

By Lemma I27T] a\ 2 > 0. The degree counting formula of (jl.ip for p £ On 
can be computed by the degree of the following specific system 



(5.20) 



A g m + a 12 p2 
A g u 2 + a u pi 
A 5 «i + piau 

+P2di2 



h 2 e u 2 



J M h 2 e^2dV g 

feie"l 
f M h ie ^dV 3 

feie"l 
J M h ie ^ldV g " 
h 2 e u 2 _ 
J M h 2 e^dV g 



o.. 
o.. 



1=0, i > 3. 



where p\ , p 2 satisfy 

(5.21) 8vr(iV - l)(pi + p 2 ) < 2a 12 p lP2 < 8vriV(pi + p 2 ). 

Easy to see (pi, p 2 , 0, .., 0) E 0at_i. Now we consider the special case of 
(|5.20p : pi = p 2 and h\ = h 2 = h. In this case, the maximum principle 
implies u\ = u 2 + c for some constant c. Since u\,u 2 are both in H l (M), 
c = 0. Then the first two equations of (|5.20p turn out to be: 

(5.22) Asu + Ol2 ^___ 1 j = o, 
where /? satisfies 

8vr(iV - 1) < a 12 p < 8ttN. 

Hence, again the Leray-Schauder degree for equation (II. ip can be reduced 
to the Leray-Schauder degree for the single equation (15.221) . By applying 
Theorem 1.2 in [13], the degree counting formulas (jl.6j) is also obtained in 
this case. This completes the proof of Theorem ll.2[ □ 
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6. Proof of Theorem 11.41 

For equation (j 1 . 3 [) . we have to show u k never blows up near the boundary 
d£l. This fact is standard, we include the argument for the convenience of 
the reader (see [34]). Since O is a bounded set with smooth boundary, there 
is a uniform constant ro such that for any point on dil, there is a ball of 
radius ro tangent to d£l at this point from the outside. Let xq £ dtt and 
B(xi,\) be a ball tangent to d£l at xq from the outside. A < ro will be 
determined later. Let 

Hi = -rrh-> ^>°> ie/ 

Jq hie Ut ax 
then the equation for u, becomes 

n 

-Am = ajjEji 1 ' , 0, iel. 

3=1 

For H{ we obviously have 

(6.1) \VlogHi(x)\<C, Viefi. 



Let y = x — x\ and 



Ui (y) = Ui(xi + A —tk), iel. 

\y\ 



Then 

-Au}(y) = YaJ^H^ + X 2 -^))^^ in A 



where f2 is the image of under the Kelvin transformation. Moreover, we 
have 

= on dn x . 



Let 



Hi(y) = ^Hi( Xl + \ 2 t^) 

\y\ \y\ 



Using (|6.ip . we see by direct computation that in a small neighborhood of 
xq, Hi is strictly decreasing in the outer normal direction to dQ x , as long 
as A is small. The smallness of the neighborhood of xq and A can both be 
represented by eo which depends on the usual constants. Thus we have the 
monotonicity of Hi in a neighborhood of the whole dQ x . Using the standard 
moving plane argument we see uf is increasing along the inner normal of 
dQ x in a small neighborhood of <9f2 A , which implies that for any sequence of 
function u k of (jl ,3[> . no blowup point for u k exists in a fixed neighborhood 
of dO,. Then the remaining part of the proof of Theorem 11.41 is the same as 
Theorem 1 1 . 1 1 and Theorem 11.21 So, the details are omitted here. □ 
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